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Abstract 

We study the Higgs sector of the Minimal Supersymmetric Standard 
Model, in the context of proton-proton collisions at LHC and SSC 
energies. We assume a relatively heavy supersymmetric particle spec- 
trum, and include recent results on one-loop radiative corrections to 
Higgs-boson masses and couplings. We begin by discussing present 
and future constraints from the LEP experiments. We then compute 
branching ratios and total widths for the neutral (h, H, A) and charged 
(i? ± ) Higgs particles. We present total cross-sections and event rates 
for the important discovery channels at the LHC and SSC. Promising 
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physics signatures are given by h — > 77, H — > 77 or or t + t~, 

A — > t + t~, and i — > followed by i7 + — > r + z^ T , which should 

allow for an almost complete coverage of the parameter space of the 
model. 



CERN-TH.6150/91 
ETH-TH/91-7 
December 1991 



1 



1 Introduction 



All available experimental data in particle physics are consistent with the 
Standard Model (SM) of strong and electroweak interactions, provided JTJ 

91 GeV < m t < 180 GeV (95% c.l.) (1) 

and 

57GeV<m^ (95% c.l.), (2) 

where m t and denote the masses of the top quark and of the SM Higgs 
boson, respectively. The lower limits on m t and are obtained from un- 
successful direct searches at the Tevatron and LEP. The upper limit on m t 
is obtained as a consistency condition of the SM, after the inclusion of ra- 
diative corrections, with the high-precision data on electroweak phenomena. 
Strong evidence for the existence of the top quark, with the quantum num- 
bers predicted by the SM, is also provided by the precise measurements of 
the weak isospin of the 6-quark. In the case of the Higgs boson, the situation 
is radically different. There is no experimental evidence yet that the mini- 
mal SM Higgs mechanism is the correct description of electroweak symmetry 
breaking. Fortunately, present and future accelerators will give decisive con- 
tributions towards the experimental solution of this problem. If the SM 
description of the Higgs mechanism is correct, LEP or the LHC and SSC 
should be able to find the SM Higgs boson and study its properties. 

Despite its remarkable successes, the SM can only be regarded as an 
effective low-energy theory, valid up to some energy scale A at which it is 
replaced by some more fundamental theory. Certainly A is less than the 
Planck scale, Mp ~ 10 19 GeV, since one needs a theory of quantum gravity 
to describe physics at these energies. However, there are also arguments, 
originating precisely from the study of the untested Higgs sector PI which 

— 1/2 

suggest that A should rather be close to the Fermi scale G F ~ 300 GeV. 
The essence of these arguments is the following. Triviality of the Xcp 4 theory, 
absence of Landau poles and perturbative unitarity imply that within the 
SM m v < 600-800 GeV. If one then tries to extend the validity of the SM 

— 1/2 

to energy scales A 3> G F , one is faced with the fact that in the SM there 
is no symmetry to justify the smallness of the Higgs mass with respect to 

2 For reviews of Higgs boson physics see, e.g., refs. |||| 
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the (physical) cut-off A. This is apparent from the fact that in the SM one- 
loop radiative corrections to the Higgs mass are quadratically divergent; it 
is known as the naturalness (or hierarchy) problem of the SM. Motivated by 
this problem, much theoretical effort has been devoted to finding descriptions 
of electroweak symmetry breaking which modify the SM at scales A ~ G F 
The likely possibility of such modifications is the reason why, when discussing 
the experimental study of electroweak symmetry breaking, one should not be 
confined to the SM Higgs, but also consider alternatives to it, which might 
have radically different signatures, and in some cases be more difficult to 
detect than the SM Higgs. Only after a thorough study of these alternatives 
can one be definite about the validity of the so-called 'no-lose theorems', 
stating that the physics signatures of electroweak symmetry breaking cannot 
be missed at LEP or the LHC and SSC. 

When considering alternatives to the minimal SM Higgs sector, it is natu- 
ral to concentrate on models which are theoretically motivated, phenomeno- 
logically acceptable and calculationally well-defined. The most attractive 
possibility satisfying these criteria is the Minimal Supersymmetric Standard 
Model (MSSM) J3|. This possibility is theoretically motivated by the fact 
that low-energy supersymmetry, effectively broken in the vicinity of the elec- 
troweak scale, is the only theoretical framework that can naturally accommo- 
date elementary Higgs bosons. The simplest and most predictive realization 
of low-energy supersymmetry is the MSSM, defined by 1) minimal gauge 
group: SU(3)c x SU(2)l X U(1)y', 2) minimal particle content: three gener- 
ations of quarks and leptons and two Higgs doublets, plus their superpartners; 
3) an exact discrete .R-parity, which guarantees (perturbative) baryon- and 
lepton- number conservation: R — +1 for SM particles and Higgs bosons, 
R — — 1 for their superpartners; 4) supersymmetry breaking parametrized 
by explicit but soft breaking terms: gaugino and scalar masses and trilinear 
scalar couplings. 

Besides the solution of the naturalness problem, there are other virtues 
of the MSSM which are not shared by many other alternatives to the SM 
Higgs and should also be recalled to further motivate our study. The MSSM 
successfully survives all the stringent phenomenological tests coming from 
precision measurements at LEP: in most of its parameter space, the MSSM 
predictions for the LEP observables are extremely close to the SM predic- 
tions, evaluated for a relatively light SM Higgs ||. This can be compared, for 
example, with the simplest technicolor models, which are ruled out by the re- 
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cent LEP data ||. Again in contrast with models of dynamical electroweak 
symmetry breaking, the MSSM has a high degree of predictivity, since all 
masses and couplings of the Higgs boson sector can be computed, at the 
tree-level, in terms of only two parameters, and radiative corrections can be 
kept under control: in particular, cross-sections and branching ratios for the 
MSSM Higgs bosons can be reliably computed in perturbation theory. Fur- 
thermore, it is intriguing that the idea of grand unification, which fails in its 
minimal non-supersymmetric implementation, can be successfully combined 
with that of low-energy supersymmetry: minimal supersymmetric grand uni- 
fication predicts a value of sin 2 Qy/ipiz) which is in good agreement with the 
measured one, and a value of the grand-unification mass which could explain 
why proton decay has escaped detection so far J7| . Finally, as a consequence 
of i?-parity, the lightest supersymmetric particle, which is typically neutral 
and weakly interacting, is absolutely stable, and thus a natural candidate for 
dark matter. 

Any consistent supersymmetric extension of the SM requires at least two 
Higgs doublets, in order to give masses to all charged quarks and leptons and 
to avoid gauge anomalies originated by the spin-1/2 higgsinos. The MSSM 
has just two complex Higgs doublets, with the following SU (3)c x SU(2)l x 
U(1)y quantum numbers (Q = T 3 l + Y): 

Hi = ~ (1,2,-1/2), H 2 = (^|) ~ (1,2, +1/2). (3) 

Other non-minimal models can be constructed, but they typically increase 
the number of parameters without correspondingly increasing the physical 
motivation. For example, the simplest non-minimal model || is constructed 
by adding a singlet Higgs field N and by requiring purely trilinear super- 
potential couplings. In this model, the Higgs sector has already two more 
parameters than in the MSSM. Folklore arguments in favour of this model are 
that it avoids the introduction of a supersymmetry-preserving mass parame- 

1/2 

ter n ~ G F and that the homogeneity properties of its superpotential recall 
the structure of some superstring effective theories. A closer look, however, 
shows that these statements should be taken with a grain of salt. First, in 
the low-energy effective theory with softly broken global supersymmetry, the 

i /2 

supersymmetric mass fi ~ G F could well be a remnant of local supersym- 
metry breaking, if the underlying supergravity theory has a suitable structure 
of interactions ||. Moreover, when embedded in a grand- unified theory, the 
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non-minimal model with a singlet Higgs field might develop dangerous insta- 
bilities [|IJ]. Also, the trilinear N 3 superpotential coupling, which is usually 
invoked to avoid a massless axion, is typically absent in string models. We 
therefore concentrate in this paper on the MSSM only. 

The previous considerations should have convinced the reader that the 
Higgs sector of the MSSM is worth a systematic study in view of the forth- 
coming hadron colliders, the LHC and SSC. To perform such a study, one has 
to deal with the rich particle spectrum of the MSSM. As discussed in more 
detail later, the Higgs sector contains one charged {H^) and three neutral 
(h, H, A) physical states. At the classical level, all Higgs boson masses and 
couplings can be expressed in terms of two parameters only, for example 
and tan/3 = vijv\. This makes the discussion more complicated than in the 
SM, where the only free parameter in the Higgs sector is the Higgs mass, 
rriip. In addition, when considering production and decay of Higgs bosons, 
the whole particle spectrum of the model has to be taken into account. As 
in the SM, the top-quark mass m t is an important parameter: barring the 
fine-tuned cases of a very light stop squark, or of charginos very close in mass 
to mz/2, the limits of eq. (1) are also valid in the MSSM |J. In contrast 
with the SM, also the supersymmetric i?-odd particles (squarks, sleptons, 
gauginos, higgsinos) can play an important role in the production and decay 



of supersymmetric Higgs bosons [[11]]. Clearly, to keep track simultaneously 
of all supersymmetric-particle masses would be a difficult (and confusing) 
task. We shall therefore concentrate, following the approach of ref. |l2 



on 



the limiting case where all supersymmetric-particle masses are heavy enough 
not to play an important role in the phenomenology of supersymmetric Higgs 
bosons. This is phenomenologically meaningful, since one can argue that a 
relatively light supersymmetric-particle spectrum is likely to give indepen- 
dent, detectable signatures at LEP or at the LHC and SSC. 

Another motivation for the present study is the recent realization [13] that 
tree-level formulae for Higgs-boson masses and couplings can receive large 
radiative corrections, dominated by the exchange of virtual top and bottom 
quarks and squarks in loop diagrams. For example, tree-level formulae would 
predict the existence of a neutral Higgs boson (h) lighter than the Z. If this 
were true, there would be a chance of testing completely the MSSM Higgs 
sector at LEP II, with no need for the LHC and SSC. However, can receive 
a large positive shift by radiative corrections, which can push h beyond the 
LEP II discovery reach. This makes the LHC and SSC important, not only 
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for a possible confirmation of a SUSY Higgs signal seen at LEP, but also for 
the exploration of the parameter space inaccessible to LEP. 

The phenomenology of the SM Higgs at the LHC and SSC 

has been intensely studied over the last years: a lot of effort was required 
to prove [p4| , |15l , at least on paper, that the combination of LEP and the 
LHC/SSC is sufficient to explore the full theoretically allowed range of SM 
Higgs masses. However, those results cannot be directly applied to the neu- 
tral states of the MSSM, since there are important differences in the cou- 
plings, and of course one needs to analyse separately the case of the charged 
Higgs. Even in the case in which all the i?-odd supersymmetric particles are 
very heavy, the Higgs sector of the MSSM represents a non-trivial extension 
of the SM case. Also several studies of the MSSM Higgs sector have already 
appeared in the literature. In particular, tree-level formulae for the MSSM 
Higgs boson masses and couplings are available, and they have already been 
used to compute cross-sections and branching ratios for representative values 
of the MSSM parameters ||. However, the existing analyses are not system- 
atic enough to allow for a definite conclusion concerning the discovery poten- 
tial of the LHC and SSC, even in the simple case of large sparticle masses. 
Also, they do not include radiative corrections to Higgs-boson masses and 
couplings. In this paper we plan to help filling these two gaps. The strategy 
for a systematic study of neutral supersymmetric Higgs bosons at the LHC 
was outlined in ref. ||12|| : however, at that time radiative corrections were 



not available, and also the 77 branching ratio was incorrectly encoded in the 
computer program. Our goal will be to see if LEP and the LHC/SSC can be 
sensitive to supersymmetric Higgs bosons in the whole (m^, tan f3) space. 

The structure of the paper is the following. In sect. 2 we review the 
theoretical structure of the Higgs sector of the MSSM, including radiatively 
corrected formulae for Higgs-boson masses and couplings. In sect. 3 we survey 
the present LEP I limits, after the inclusion of radiative corrections, and the 
plausible sensitivity of LEP II. In sect. 4 we present branching ratios and 
widths of neutral and charged supersymmetric Higgs bosons. In sect. 5 we 
compute the relevant cross-sections at the LHC and SSC, and in sect. 6 we 
examine in some detail the most promising signals for discovery. Finally, 
sect. 7 contains a concluding discussion of our results and of prospects for 
further work. 



6 



2 Higgs masses and couplings in the MSSM 



For a discussion of Higgs-boson masses and couplings in the MSSM, the 
obvious starting point is the tree-level Higgs potential [|J 

m\ \H^ 2 + m\ \H 2 \ 2 + m\ {H X H 2 + h.c.) 



Vn 



+\g 2 (H\aH 2 + HlaH,) 2 + \g' 2 (\H 2 \ 2 - \H X \ 



(4) 



where m^,m|,m| are essentially arbitrary mass parameters, g and g' are 
the SU{2) and U(l) coupling constants, respectively, and a are the Pauli 
matrices. SU (2) indices are left implicit and contracted in the obvious way. 
It is not restrictive to choose m\ real and negative, and then the vacuum 
expectation values v\ = (H®) and v 2 = (H 2 ) real and positive. 

The physical states of the MSSM Higgs sector are three neutral bosons 
(two CP-even, h and H, and one CP-odd, A) and a charged boson, H ± . 
A physical constraint comes from the fact that the combination (vf + v 2 ), 
which determines the W and Z boson masses, must reproduce their mea- 
sured values. Once this constraint is imposed, in the Born approximation 
the MSSM Higgs sector contains only two independent parameters. A con- 
venient choice, which will be adopted throughout this paper, is to take as 
independent parameters m^, the physical mass of the CP-odd neutral bo- 
son, and tan/5 = v 2 /vi, where v x gives mass to charged leptons and quarks 
of charge —1/3, v 2 gives mass to quarks of charge 2/3. The parameter 
is essentially unconstrained, even if naturalness arguments suggest that it 
should be smaller than 0(500 GeV), whereas for tan/? the range permitted 
by model calculations is 1 < tan/3 ^ 

At the classical level, the mass matrix of neutral CP-even Higgs bosons 
reads 



m° r y 



cot (3 
-1 



-1 

tan f3 



m A z /tan/3 -1 
~2~ + V -1 cot/3 



m t 



and the charged-Higgs mass is given by 



2 2,2 

m H ± = m w + m A . 



sin 2/3 (5) 



(6) 



From eq. @, one obtains 



m 



h,H 



2 i 2 _ 

m A + m z T 



77r 



+ m|) 2 — Am A m 2 z cos 2 2/3 



(7) 
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dd, ss, bb 

e + e~ , , t + t~ 


uu, cc, tt 




h 


— sin a/ cos (3 


cos a/ sin (3 


sin (/3 — a) 


H 


cos a / cos (3 


sin a/ sin /5 


cos (/3 — a) 


A 


— «7 5 tan (3 


—275 COt /3 






Table 1: Correction factors for the couplings of the MSSM neutral Higgs 
bosons to fermion and vector boson pairs. 



and also celebrated inequalities such as m w ,m A < m H ±, m h < m z < run, 
rrih < rriA < ran- Similarly, one can easily compute all the Higgs-boson 
couplings by observing that the mixing angle a, required to diagonalize the 
mass matrix (^j), is given by 

2 2 ,_ 

cos la = — cos 1(3 —75————^ , -- < a < 0. (8) 
mjj — m h 2 

For example, the couplings of the three neutral Higgs bosons are easily ob- 
tained from the SM Higgs couplings if one multiplies them by the a- and 
/3-dependent factors summarized in table 1. The remaining tree- level Higgs- 
boson couplings in the MSSM can be easily computed and are summarized, 
for example, in ref. |J. An important consequence of the structure of the 
classical Higgs potential of eq. is the existence of at least one neutral 
CP-even Higgs boson, weighing less than or about mz and with approxi- 
mately standard couplings to the Z. This raised the hope that the crucial 
experiment on the MSSM Higgs sector could be entirely performed at LEP 
II (with sufficient centre-of-mass energy, luminosity and 6-tagging efficiency), 
and took some interest away from the large hadron collider environment. 



8 



However, it was recently pointed out [[HJ that the masses of the Higgs 
bosons in the MSSM are subject to large radiative corrections, associated 
with the top quark and its SU{2) and supersymmetric partners^. Several 
papers p0|-|2"3|1 have subsequently investigated various aspects of these cor- 
rections and their implications for experimental searches at LEP. In the rest 
of this section, we shall summarize and illustrate the main effects of radiative 
corrections on Higgs-boson parameters. 

As far as Higgs-boson masses and self-couplings are concerned, a conve- 
nient approximate way of parametrizing one-loop radiative corrections is to 
substitute the tree-level Higgs potential of eq. (f|) with the one-loop effective 
potential, and to identify Higgs-boson masses and self-couplings with the 
appropriate combinations of derivatives of the effective potential, evaluated 
at the minimum. The comparison with explicit diagrammatic calculations 
shows that the effective potential approximation is more than adequate for 
our purposes. Also, inspection shows that the most important corrections 
are due to loops of top and bottom quarks and squarks. At the minimum 
(Hf) = vi, (H®) = t> 2 , (Hi ) = (H2) = 0, and neglecting intergenerational 
mixing, one obtains for the top and bottom quark and squark masses the 
familiar expressions 

2 l2 2 2 7,2 2 /n\ 

m t =h t v 2 , m h = h b v 11 (9) 

2 2,1/2, 2 \ 1 1 _2 no 

m h,2 = m t + 2 y m Q + m u) + cos 2(3 
± 



~(m 2 Q - ml) + ^(8m^ - 5m|) cos 2(3 



+ m| (A t + /i cot (3) 



2 1 / 2 1 2 \ 

m b + ^( m Q + m D) 



-m 



I cos 2(3 



(m 2 Q - m 2 D ) 



~ j^( Am w cos 2/3 



+ m 2 (Ab + fj, tan (3)' 



3 Previous studies j[9[ either neglected the case of a heavy top quark, or concentrated 
on the violations of the neutral-Higgs mass sum rule without computing corrections to 
individual Higgs masses. 
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In eqs. (^) to (|TT]), h t and /i& are the top and bottom Yukawa couplings, 



and mQ,mu,m£i are soft supersymmetry-breaking squark masses. The pa- 
rameters At, A b and //, which determine the amount of mixing in the stop 
and sbottom mass matrices, are defined by the trilinear potential terms 
h t A t (t L t c L H$ - ~b L t c L H+) + h.c, h h A h (b L ¥ L R\ - t L b c L H{) + h.c. and by 
the superpotential mass term ^(H^H® — H^H^), respectively. 

To simplify the discussion, in the following we will take a universal soft 
supersymmetry-breaking squark mass, 

2 2 2 2 /io\ 

m Q = m U = m D = m qi ( 12 ) 

and we will assume negligible mixing in the stop and sbottom mass matrices, 

A t = A b = /2 = 0. (13) 
Formulae valid for arbitrary values of the parameters can be found in refs. 



2^,^], but the qualitative features corresponding to the parameter choice 
of eqs. (|12|) and (|13|) are representative of a very large region of parameter 
space. In the case under consideration, and neglecting D-term contributions 
to the field-dependent stop and sbottom masses, the neutral CP-even mass 
matrix is modified at one loop as follows 

Ml=(M%)\(^ ° 2 ), (14) 

where 



A 2 



3g 2 m 4 b 



2 1 

m~ m~ 



lo7r z 77% cos z p ml 

a9 3q 2 m 4 , m 2 -m 2 r 

167r 2 77% sin p vrq 

From the above expressions one can easily derive the one-loop-corrected 
eigenvalues and m#, as well as the mixing angle a associated with the 
one-loop-corrected mass matrix ([14]). The one-loop-corrected charged Higgs 



mass is given instead by 

™%± =m 2 w + m 2 A + A 2 , (17) 
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where, including D-term contributions to stop and sbottom masses, 



A 2 



3£ 2 

647r 2 sin 2 (3 cos 2 /3to| 



w 



(ml — m^r cos 2 /3) (to 2 — mfy sin 2 (3) - 



m z 



mi 



f(ml)-f(ml) 



+ 



2 2 

m t ml 

2 2~~ 

ml — m~ 

*2 02 



/(to|)-/(to?J 



2m 2 m 2 



mr 



kg) 



and 

log— -lj. (19) 

The most striking fact in eqs. (|I4D— (|I9|) is that the correction A 2 , is propor- 
tional to (mj/myy)- This implies that, for m t in the range of eq. (HD, the 
tree-level predictions for m h and m H can be badly violated, and so for the 
related inequalities. The other free parameter is m^, but the dependence on 
it is much milder. To illustrate the impact of these results, we display in fig. 1 
contours of the maximum allowed value of m^ (reached for m& — > oo), in the 
(mt, tan/3) and (m t ,mg) planes, fixing m§ = 1 TeV and tan/3 = m t /m , re- 
spectively. In the following, when making numerical examples we shall always 
choose the representative value mq = 1 TeV. To plot different quantities of 
physical interest in the (to^, tan/5) plane, which is going to be the stage of 
the following phenomenological discussion, one needs to fix also the value of 
m t . In this paper, whenever an illustration of the m t dependence is needed, 
we work with the two representative values m t = 120, 160 GeV, which are 
significantly different but well within the range of eq. ([I]). Otherwise, we 
work with the single representative value m t = 140 GeV. As an example, we 
show in figs. 2-4 contours of constant m^, mjj, and to#± in the (m^, tan/3) 
plane. Here and in the following we vary m^ and tan (3 in the ranges 

< m A < 500 GeV, 1 < tan/3 < 50 . (20) 

The effective-potential method allows us to compute also the leading cor- 
rections to the trilinear and quadrilinear Higgs self-couplings. A detailed 
discussion and the full diagrammatic calculation will be given elsewhere. 



11 



Here we just give the form of the leading radiative corrections to the trilin- 
ear hAA, HAA, and Hhh couplings, which will play an important role in 
the subsequent discussion of Higgs-boson branching ratios. One finds |23],E1 



^hAA — tfiAA + ^^hAA > ^HAA — ^HAA + ^^HAA , ^Hhh — ^°Hhh + ^^Hhh , 

(21) 

where 

Kaa = -T^h cos2/3sin(/3 + a) , (22) 
z cos t>w 

*haa= ^^cos2/5cos(/? + a), (23) 

Z COS u\y 

^Hhh = - n ^ m f [2 sm((3 + a) sin 2a - cos(fi + a) cos 2a] , (24) 
2 cos vw 

and, neglecting the bottom Yukawa coupling and the .D-term contributions 
to squark masses 

igm z 3g 2 cos 2 9 W cos a cos 2 p m\ m\ + m 2 
AA /tA A = —r- — -„ — 2 —3-5 — — — log — —J , (25) 

Z COS vyy oTT Sill p TTl w TTl t 

igm z 3g 2 cos 2 6 W sin a cos 2 /3 m\ m\ + m 2 

AA^ AA = — r^— — log , 26 

2 costly on 2 sm p 77% mf 

igm z 3g 2 cos 2 9 w cos 2 a sin a m\ ( m 2 - + m 2 m~ 
AA HWl = -77—2 ^-2 —r3 — — — 31 °S — —2 2 Z2j^2 

Z COS U\y o7T Sin p ^^(/ \ TTlf 771= + 771^ 

(27) 

Notice that, besides the obvious explicit dependence, in eqs. (pl])-([27|) there is 
also an important implicit dependence on m t and rriq, via the angle a, which 
is determined from the mass matrix of eqs. (|T^) — (|TBD • We also emphasize that 
neglecting the D-terms in the stop and sbottom mass matrices is guaranteed 
to give accurate results only for m t ^> mz- For m t ~ m z , one should make 
sure that the inclusion of Z)-terms does not produce significant modifications 
of our results. In the case of the h and H masses, and of the mixing angle 
a, complete formulae are available, and this check can be easily performed. 
In the case of the hAA, HAA and Hhh couplings, complete formulae are 
not yet available. For the phenomenologically most important coupling at 
the LHC and SSC, Xnhh, we have explicitly checked that the inclusion of 



D-terms does not produce important modifications of our results. 
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Finally, one should consider Higgs couplings to vector bosons and fermions. 
Tree-level couplings to vector bosons are expressed in terms of gauge cou- 
plings and of the angles j3 and a. The most important part of the radiative 
corrections is taken into account by using one-loop-corrected formulae to de- 
termine a from the input parameters. Other corrections are at most of order 
g 2 mf/m'^ v and can be safely neglected for our purposes. Tree- level couplings 
to fermions are expressed in terms of the fermion masses and of the angles 
/3 and a. In this case, the leading radiative corrections can be taken into ac- 
count by using the one-loop-corrected expression for a and running fermion 
masses, evaluated at the scale Q which characterizes the process under con- 
sideration. This brings us to the discussion of the renormalization group 
evolution of the top and bottom Yukawa couplings in the MSSM. As bound- 
ary conditions, we assume as usual that m t (m t ) = m t and m&(m{,) = nif,, with 
rrib = 4.8 GeV and m t numerical input parameters. As stated in the Intro- 
duction, we assume in this paper that all supersymmetric particles are heavy. 
Then, since we want to compute Higgs-boson production cross-sections and 
branching ratios, we are interested in the standard renormalization group 
evolution of h t {Q) [hb(Q)) from Q = m t [Q = nib] to Q — m H ±,rriH, which 
is dominated by gluon loops. 

To illustrate the behaviour of the Higgs couplings to vector bosons and 
fermions, as functions of the input parameters, we show in figs. 5-7 contours 
in the (m^, tan (5) plane of some of the correction factors appearing in table |I[ 



3 LEP limits and implications 

In this section, we briefly summarize the implications of the previous results 
on MSSM Higgs boson searches at LEP I and LEP II. Partial results were 
already presented in refs. |21,22| . 

As already clear from tree-level analyses, the relevant processes for MSSM 
Higgs boson searches at LEP I are Z — > hZ* and Z — > hA, which play 
a complementary role, since their rates are proportional to sin 2 (/3 — a) and 
cos 2 (/5— a), respectively. An important effect of radiative corrections |23j is to 
allow, for some values of the parameters, the decay h — > AA, which would be 
kinematically forbidden according to tree-level formulae. Experimental limits 
which take radiative corrections into account have by now been obtained by 



the four LEP collaborations using different methods to present and 
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analyse the data, and different ranges of parameters in the evaluation of 
radiative corrections. A schematic representation of the presently excluded 
region of the (m^, tan (3) plane, for the standard parameter choices discussed 
in sect. 2, is given in fig. 8, where the solid lines correspond to our naiveF] 



extrapolation of the exclusion contour given in the first of refs. ||25|| . For 
a discussion of the precise experimental bounds, we refer the reader to the 
above-mentioned experimental publications. 

The situation in which the impact of radiative corrections is most dra- 
matic is the search for MSSM Higgs bosons at LEP II. At the time when 
only tree- level formulae were available, there was hope that LEP could com- 
pletely test the MSSM Higgs sector. According to tree-level formulae, in fact, 
there should always be a CP-even Higgs boson with mass smaller than (h) 
or very close to (H) mz, and significantly coupled to the Z boson. However, 
as should be clear from the previous section, this result can be completely 
upset by radiative corrections. A detailed evaluation of the LEP II discov- 
ery potential can be made only if crucial theoretical parameters, such as the 
top-quark mass and the various soft supersymmetry-breaking masses, and ex- 
perimental parameters, such as the centre-of-mass energy, the luminosity and 
the 6-tagging efficiency, are specified. Taking for example \fs = 190 GeV, 
m t ^ 110 GeV, and our standard values for the soft supersymmetry-breaking 
parameters, in the region of tan/3 significantly greater than 1, the associated 
production of a Z and a CP-even Higgs can be pushed beyond the kine- 
matical limit. Associated hA production could be a useful complementary 
signal, but obviously only for rrih + vtla < y/s. Associated HA production 
is typically negligible at these energies. To give a measure of the LEP II 
sensitivity, we plot in fig. 8 contours associated with two benchmark val- 
ues of the total cross-section a(e + e~ — > hZ*,HZ*,hA,HA). The dashed 
lines correspond to a = 0.2 pb at y/s = 175 GeV, which could be seen as 
a rather conservative estimate of the LEP II sensitivity. The dash-dotted 
lines correspond to a = 0.05 pb at y/s = 190 GeV, which could be seen as 
a rather optimistic estimate of the LEP II sensitivity. In computing these 
cross-sections, we have taken into account the finite Z width, but we have 



4 We fitted the experimental exclusion contours, corresponding to m f = 140 GeV 
and the other parameters as chosen here, with two numerical values for T(Z — > hZ*) 
and Y(Z — > hA). We have then computed radiative corrections for the two values of mt 
considered here, assuming that the variations in experimental efficiencies are small enough 
not to affect our results significantly. 
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neglected initial state radiation, which leads to suppression near threshold. 
A more accurate estimate of the LEP II sensitivity can be found in ref. |26 . 
Of course, one should keep in mind that there is, at least in principle, the 
possibility of further extending the maximum LEP energy up to values as 
high as a/s ~ 230-240 GeV, at the price of introducing more (and more 
performing) superconducting cavities into the LEP tunnel [p7f| . 

In summary, a significant region of the parameter space for MSSM Higgses 
could be beyond the reach of LEP II, at least if one sticks to the reference 
centre-of-mass energy y/s ^190 GeV. The precise knowledge of this region 
is certainly important for assessing the combined discovery potential of LEP 
and LHC/SSC, but it does not affect the motivations and the techniques of 
our study, devoted to LHC and SSC searches. Whether or not a Higgs boson 
will be found at LEP in the future, we want to investigate the possibilities of 
searching for all the Higgs states of the MSSM at large hadron colliders, in 
the whole region of parameter space which is not already excluded at present. 
Even if a neutral Higgs boson is found at LEP, with properties compatible 
with the SM Higgs boson within the experimental errors, it will be impossible 
to exclude that it belongs to the MSSM sector. The LHC and SSC could then 
play a role in investigating its properties and in looking for the remaining 
states of the MSSM. 

Similar considerations can be made for charged-Higgs searches at LEP II 
with t/s ^ 190 GeV. In view of the /3 3 threshold factor in a(e + e~ — > H + H~), 
and of the large background from e + e~ — *> W + W~, it will be difficult to 
find the H ± at LEP II unless m#± ^ m^y, and certainly impossible unless 
m H± < \/s/2. We also know |23|,[2(J that for generic values of the parameters 
there are no large negative radiative corrections to the charged-Higgs mass 
formula, eq. ([]). A comparison of figs. 4 and 8 indicates that there is very 
little hope of finding the charged Higgs boson of the MSSM at LEP II (or, 
stated differently, the discovery of a charged Higgs boson at LEP II would 
most probably rule out the MSSM). 



4 Branching ratios 
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4.1 Neutral Higgs bosons 

The branching ratios of the neutral Higgs bosons of the MSSM were sys- 
tematically studied in ref. 0], using the tree-level formulae for masses and 
couplings available at that timeQ (previous work on the subject is summa- 
rized in ref. @). Here we present a systematic study which includes the 
radiative corrections described in sect. 2. As usually done for the SM Higgs 
boson, we consider the two-body decay channels 

h,H,A — >cc, bb, ti,T + T~,gg } ^ } W*W* } Z*Z*, Z<y . (28) 

For consistency, we must also consider decays with one or two Higgs bosons 
in the final state 

h — > AA, H -> hh, AA, ZA, A -> Zh. (29) 

On the other hand, we neglect here possible decays of MSSM Higgs bosons 
into supersymmetric particles: as previously stated, we consistently assume 
a heavy spectrum of i?-odd particles, so that only .R-even ones can be kine- 
matically accessible in the decays of h, H, and A. We perform our study in 
the framework of MSSM parameter space, with the representative parameter 
choices illustrated in sect. 2. The effects of changing the mass of the top 
quark, and the sensitivity to squark masses in the high-mass region, will also 
be briefly discussed. 

The partial widths for the decays of eq. that correspond to tree- 

level diagrams can be obtained from the corresponding formulae for the SM 
Higgs boson (for a summary, see ref. ||), by simply multiplying the vari- 
ous amplitudes by the supersymmetric correction factors listed in table 1. 
For decays that are described by loop diagrams, however, in the MSSM 
one has to include some contributions that are absent in the SM. Diagrams 
corresponding to the exchange of R-odd supersymmetric particles give negli- 
gible contributions to the corresponding partial widths, in the limit of heavy 
supersymmetric-particle masses that we have chosen for our analysis (in ac- 
cordance with intuitive ideas about decoupling). One must also include the 
charged-Higgs loop contributions to the 77 and Z7 final states. When con- 
sidering instead the processes of eq. (|2"9"D, we improve the tree-level formulae 

5 Also, the partial widths for the decays h, H, A — > 77 were affected by numerical errors. 
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of ref. not only with the self-energy corrections to the mixing angle a, 
but also with the vertex corrections of eqs. (|2T|) -(|2"7l ) . 

QCD [^] and electroweak [2~9~| radiative corrections to the fermion-anti- 



fermion and the WW, ZZ channels have been recently computed for the 
SM Higgs boson, <p. They have been found to be small (less than ~ 20%), 
with the exception of the QCD corrections to the decays into charm and 
bottom quark pairs, which are large because of running-quark-mass effects. 
We then included the QCD corrections as described in ref. [|14[]. One may 



also wonder whether running-mass effects induced by the large top Yukawa 
coupling could give further important effects. However, one can easily see 
that these effects give corrections which are certainly less than 20%. 

The QCD correction to tp — > 77 is also available, and known to be negli- 
gibly small Sizeable QCD corrections are found, however, for the decay 



tp — > gg [51]. Although this effect is not important for the branching ratio 
study, since ip ^ gg is neither the dominant decay mode nor a useful channel 
for detection, it still has to be included in the production cross-section of h 
via the two-gluon fusion mechanism. 

Another general and well-known property of the MSSM is that the self- 
interactions of the Higgs bosons are controlled, modulo the logarithmic cor- 
rections discussed in sect. 2, by the SU{2) and U{1) gauge couplings. There- 
fore, the total widths of all MSSM Higgs bosons, displayed in fig. 9, stay 
below 10 GeV in the whole parameter space we have considered. 

The most important branching ratios for the neutral MSSM Higgs bosons 
are shown, as a function of the mass of the decaying particle, in figs. 10-12. 
To avoid excessive proliferation of figures, we consider the two representative 
values 

tan p = 1.5, 30, (30) 

and for each of these we vary tua between the experimental lower bound of 
fig. 8 (m A ~ 59 GeV for tan/? = 1.5, m A ~ 44 GeV for tan (3 = 30) and the 
upper bound of eq. (^3), assuming m t = 140 GeV and nig = 1 TeV. 

We consider first the branching ratios of h (fig. 10). We can clearly see 
the effect of radiative corrections on the allowed range of rrih for the given 
values of tan (3. For itla ^ 25 GeV, the decay h — > AA can be kinematically 
allowed and even become the dominant mode. This decay channel was im- 
portant at LEP I, but since the corresponding region of parameter space is 
already excluded by experiment, this decay mode does not appear in fig. 10. 
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The dominant decay mode is then h — > 66, whereas the t + t~ mode has a 
branching ratio of about 8% throughout the relevant part of the parameter 
space. In fig. 10, one immediately notices the rather steep slopes for the 
cc and 77 branching ratios plotted versus m^, with larger effects for larger 
values of tan/3: their origin can be understood by looking at figs. 2 and 5-7, 
which show how rrih and the h couplings to heavy fermions and vector bosons 
vary in the (m^,tan/3) plane. 

If the SM Higgs boson is in the intermediate mass region, m v = 70-140 
GeV, at large hadron colliders a measurable signal can be obtained via the 77 
mode. Since the mass of the light Higgs h is indeed below or inside this region, 
the 77 mode is also crucial for the MSSM Higgs search. Furthermore, the 
77 branching ratio as a function of the Higgs mass exhibits a rather peculiar 
behaviour, not only for h but also for H and A, so a more detailed discussion 
is in order. The partial width is given by 



p/ JL \ Q V m l 

r( ^ 77) = T024^^T 
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where <fi = h,H,A and i — f, W, H, f, \ indicates the contributions from 
ordinary fermions, charged gauge bosons, charged Higgs bosons, sfermions 
and charginos, respectively. The functions I\(t1) are given by 

I* = Ff /2 {rj)N cf e) B$ , 
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It = F„(rpN^Rf^, 

It = Ft /2 (r*)4^ (32) 

where N c f is 1 for (s)leptons and 3 for (s)quarks, and the subscripts of the 
complex functions Ff, 2 (r), F[j 2 (t), F (t), and Fi(r), which were calculated 
in ref. [[32]] , indicate the spin of the particles running in the loop. In the case 
of spin-1/2 particles, the contribution is different for CP-even and CP-odd 
neutral Higgses. The symbols Ftf denote the appropriate correction factors 
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for the MSSM Higgs couplings: for i = f, W they are given in table 1, for 
i = H, f,x they can be found, for example, in Appendix C of ref. |J. The W 
contribution dominates the h — > 77 decay rate. The function Fi is large at 
and above r = 1. For the W contribution r = 4m^/m^ > 1, and increasing 
m h gives increasing values of F 1 . The steep dependence of the branching ratio 
on mh is a consequence of the fast change of sin 2 (/3 — a) as is increased 
for fixed tan/?. This is further enhanced by the fact that the large interval 
100 GeV ^ rriA < 500 GeV is mapped into a very small interval (a few GeV) 
in m/j. We elucidate this effect by plotting in fig. 13 the branching ratios of 
h as a function of m^, for the same values of the parameters as in fig. 10. 
We can see that the tip of the gg, cc and 77 curves in fig. 10 is mapped into 
a long plateau in fig. 13. We can also observe that in a large region of the 
parameter space the h — > 77 branching ratio has a value somewhat smaller 
than (but comparable to) the corresponding branching ratio for a SM Higgs 
of mass m/j. This is due to the fact that all the h couplings tend to the SM 
Higgs couplings for rriA mz\ however, for the h couplings to fermions the 
approach to the asymptotic value is much slower than for the h couplings 
to vector bosons, as can be seen from figs. 5 to 7. In fig. 13, the branching 
ratios for the W*W* and Z*Z* decays are also plotted, whereas they were 
omitted in fig. 10 in order to avoid excessive crowding of curves. However, 
for our parameter choice they have little interest at large hadron colliders, 
because of the small production rates and the large backgrounds. 

The branching ratios of the heavy Higgs boson H, depicted in fig. 11, 
have a rather complicated structure. We make here four remarks. 

i) The 77 mode has a steeply decreasing branching ratio with increasing 
iriHi except at small values of tan/3 and at the lower kinematical limit of 
rriH, where one or more of the A A, ZA and hh decay channels are open. 
The steep fall of the 77 branching ratio at large values of tan/? can be 
easily understood. The partial width T(H — > 77) is dominated by the W 
contribution, proportional to cos 2 (/3 — a). As we can see in fig. 7, cos 2 (/3 — 
a) decreases very fast, for increasing m#, at fixed values of tan/?. This 
steep decrease is slightly compensated by the increase of Fi(tw) at Tw < 1, 
which has a peak at the W threshold m# = 2m w Another peak in the 77 
branching ratio is obtained, for small values of tan/3, at itlh = 2m t , where 
the top-quark loop gives the dominant contribution. 

ii) The complicated structure in the H branching ratio curves is mainly 
due to the H —>■ hh channel. For inn < "2rn t , and not too high values 
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of tan/3, this decay mode is dominant whenever kinematics allows. This 
channel is always open at the lower kinematical limit of Mh- Increasing Mh 
a little bit, however, it may become strongly suppressed, because for small 
increasing values of wia the value of rah rises faster than that of m#, so that 
the channel can become kinematically closed. Obviously, for sufficiently high 
values of m# the channel is always open. At high values of tan/3, the mass 
region at the lower kinematical limit where H — ► hh is open becomes smaller 
and smaller, explaining the presence of the almost vertical line in fig. 11. A 
further structure is present in this decay channel due to the coupling factor 
^Hhh [see eqs. fl2"3|) and fl2"T|)1. There are relatively small values of ran at 
which \uhh accidentally vanishes. Furthermore, for very large values of ran 
and tan/3 one has a ~ 0, (3 ~ ir/2, and therefore Xhuh — 0- Unfortunately, 
even when it is dominant, this mode has very large backgrounds, so it seems 
unlikely to give a measurable signal at large hadron colliders. The H — > AA 
mode is kinematically allowed only for values of ra^ below 50-60 GeV, in 
which case it can have a large branching ratio, competing with the one for 
H — ► hh. The H — > ZA mode is kinematically allowed only in the region of 
parameter space which is already excluded by the LEP I data. 

iii) H can decay at tree level into ZZ — > l + l~l + l~, which is the 'gold- 
plated' signature for the SM Higgs boson. Unfortunately, in the case of H 
the branching ratio is smaller, and it decreases fast with increasing tan/3 
and/or m^. For small tan/3 and Im^ < ran < 2rat, this mode is suppressed 
by the competition with H — > hh, and this effect is further enhanced by 
the inclusion of the radiative correction of eq. p7|), which typically gives an 
additional 50% suppression. Nevertheless, as we shall see in the next section, 
the four-lepton channel can give a measurable signal in some small region of 
the parameter space. 

iv) The decay into ti is dominant above threshold at moderate values of 
tan/3. But above tan/5 ~ 8 or so the bb mode remains dominant and t + t~ 
has the typical ~ 10% branching ratio. 

Finally, we discuss the branching ratios of A, shown in fig. 12. The 77 one 
is always small, although at small tan /3 and slightly below the top threshold, 
raA ~ 2m t , it reaches a value ~ 8 x 10~ 4 , which may give a measurable signal 
in a small island of the parameter space. The behaviour of the A — > 77 
branching ratio can be easily understood by taking into account that the 
partial width is dominated by the top loop contribution. Two features are 
important here. First, the function -Fi/2(t) appearing in eq. fl32|) has a strong 
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enhancement at r ~ 1. Furthermore, the ttA coupling gives a suppression 
factor l/(tan/3) 2 for increasing values of tan/3. 

At smaller values of and tan f3, there is a substantial branching ratio 
to Zh, which however does not look particularly promising for detection at 
large hadron colliders, because of the very large Zbb background. We can see 
that all the dominant decay modes of the A boson correspond to channels 
which are overwhelmed by very large background, except perhaps the t + t~ 
mode, which, as we shall see in the next section, may give a detectable signal 
for very high values of tan f3. 

We have studied the neutral Higgs branching ratios also at m t = 120, 160, 180 GeV. 
Increasing the top mass has two major effects. First, the maximum value of 
m/j increases (see fig. 1). Next, owing to the increased value of the top 
threshold, the structure generated by the opening of the top decay channel 
is shifted to higher mass values. We also note that varying rriq in the range 
0.5-2 TeV has negligible effects on the branching ratio curves of figs. 10-12. 
Finally, if one chooses m t and rriq so large that rrih > 130 GeV, the WW* 
and Z*Z* branching ratios can become relevant also for h. 

4.2 Charged Higgs boson 

In the case of the charged Higgs boson, we considered only the two-body 
decay channels 

H + -> cs, r+z/ r , tb, W + h. (33) 

Tree-level formulae for the corresponding decay rates can be found, for in- 
stance, in ref. ||. Loop- induced decays such as H + — > W +, y, W + Z have very 
small branching ratios and are not relevant for experimental searches at 



the LHC and SSC Radiative corrections to the charged-Higgs-boson mass 



formula were included according to eqs. (|T7|) and ([18|) . The H + W~h cou- 
pling, proportional to cos(/3 — a), was evaluated with the one-loop corrected 
value of a. The leading QCD corrections to the H + bt and H + sc vertices were 
parametrized, following refs. |34j], by running quark masses evaluated at a 
scale Q ~ rrin± . The resulting branching ratios for the charged Higgs boson 
are displayed in fig. 14, for tan f3 = 1.5, 30 and the standard parameter choice 
m t = 140 GeV, rriq = 1 TeV. One can see that the dominant factor affecting 
the branching ratios is the = m t + rrib threshold. Above threshold, the 
tb mode is dominant for any value of tan (3 within the bounds of eq. (|20| ) . 
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Below threshold, the dominant mode is t + u t , with the competing mode cs 
becoming more suppressed for higher values of tan/3. For small values of 
tan/3, the W + h decay mode can also be important, and even dominate, in 
a limited m H ± interval, if the W + h threshold opens up before the tb one. 
The exact position of the two thresholds on the m#± axis depends of course 
on tan/5, m t , and rriq. It is just a numerical coincidence that in fig. 14a 
the two thresholds correspond almost exactly. For increasing values of m#± 
and tan/?, the numerical relevance of the W + h branching ratio rapidly dis- 
appears, because of the cos 2 (/3 — a) suppression factor in the corresponding 
partial width. 

The total charged Higgs boson width is shown, as a function of m#± and 
for tan/3 = 1.5, 3, 10, 30, in fig. 9d. Again one can see the effects of the tb 
threshold, and also the tan /^-dependence of the couplings to fermions. In 
any case, the charged Higgs width remains smaller than 1 GeV for m#± < 
m t + m&, and smaller than 10 GeV for m H ± < 500 GeV. 



5 Neutral-Higgs production cross-sections 

There is only a limited number of parton-level processes which can give in- 
teresting rates for the production of the MSSM neutral Higgs bosons (0 = 
h, H, A) at proton-proton supercolliders: 

9 + 9 -> , (34) 

q + q-^q + q + W* + W*-^q + q + <p, (35) 

g + g or q + q — > b + b + , (36) 

g + g or q + q —>• t + t + (j) , (37) 

q + q^W(Z) + <j), (38) 

where q denotes any quark flavour. These processes are controlled by the 
Higgs couplings to heavy quarks and gauge bosons, whose essential fea- 
tures were summarized in table 1. We briefly discuss here the corresponding 
cross-sections and the status of their theoretical description, emphasizing the 
features which are different from the SM case. We shall always adopt the 
HMRSB structure functions [0511 with = 190 MeV. 
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Gluon fusion. In the SM, gg — > <fi [36] is the dominant production 
mechanism, the most important diagram being the one associated with the 
top-quark loop. In the MSSM, this is not always the case, since the correction 
factors of table 1 give in general suppression for the top contribution and 
enhancement for the bottom one, and stop and sbottom loops could also 
play a role. 

The leading-order amplitudes for the gluon-fusion processes are deter- 
mined by the functions of eqs. (|3~2"D, with top, bottom, stop and sbottom 
intermediate states. For rriq = 1 TeV, the squark contributions are very 
small, owing to the suppression factor m|/m~ in the corresponding Iq func- 
tions. For large values of tan /3, the bottom contribution can compete with 
the top one and even become dominant. 

QCD corrections to the gluon-fusion cross-section were recently evaluated 
in ref. PTfl , for a SM Higgs in the mass region below the heavy-quark thresh- 



old. In this region, QCD corrections increase the top contribution by about 
50%. To a good approximation, the bulk of QCD corrections can be taken 
into account by performing the replacement 



1 + 1- + ^ as 
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(39) 



at the renormalization scale Q = m^. This calculation, unfortunately, is 
not valid above the heavy-quark threshold, a region which is relevant for 
the bottom contribution and for the top contribution to H, A production, 
when rriH,mA > 2mt- Even below the heavy-quark threshold, the SM QCD 
corrections are applicable to h and H production, but not to A production, 
because of the additional 75 factor appearing at the Aqq vertex. In view of 
this not completely satisfactory status of QCD corrections, we calculate, con- 
servatively, the top contribution without QCD corrections. However, when 
discussing the detectability of the different physics signals, we shall take into 
account the results of ref. pi |, when applicable. In the case of the bottom 



contribution, we use the running m&, which leads to suppression. 

In fig. 15 we display cross-sections for gg — > (0 = h, H, A), as functions 
of m^, for tan (3 = 1.5, 3, 10, 30 and for LHC and SSC energies. The SM Higgs 
cross-section is also shown for comparison. For large values of tan (3 and not 
too high values of m^, the cross-sections can be enhanced with respect to the 
SM value. This effect is due to the enhanced bottom-quark contribution, as 
apparent from table 1 and fig. 5. The fast disappearance of this effect for 
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increasing Higgs masses is due to the fast decrease of the function Fy 2 (jf) 

as — > 0. When the neutral Higgs couplings to fermions are SM-like, the 
gluon-fusion cross-sections approach the SM value, and are always dominated 
by the top contribution. The changes in the slopes of the curves in fig. 15 are 
due to the competing top and bottom contributions. In particular, one can 
notice an important threshold effect, for ~ 2m t , in the process gg — > A, 
which can bring the corresponding cross-section above the SM value for low 
tan (5. 

As a final remark, we notice that the LHC and SSC curves in fig. 15 have 
very similar shapes, with a scaling factor which is determined by the gluon 
luminosity and uniformly increases from ~ 2.5 at ~ 100 GeV to ~ 5 at 

~ 500 GeV. 

W fusion. In the SM case, the VT-fusion mechanism J37j can compete 
with the gluon-fusion mechanism only for a very heavy (m^ 500 GeV) 
Higgs boson, owing to the enhanced WlWl^P coupling and to the relative 
increase of the quark number-densities. In the MSSM, the correction factors 
for the couplings to vector boson pairs (see table 1 and fig. 7) are always 
smaller than 1, so that the MSSM IV-fusion cross-sections are always smaller 
than the SM one. 

We illustrate this in fig. 16, where iy-fusion cross-sections for h and H 
are displayed, for the same tan/3 and y^i values as in fig. 15. For both 
h and H, the SM cross-section is approached from below in the regions of 
parameter space where sm 2 (/3 — a) — > 1 and sin 2 (/3 — a) — > 0, respectively. In 
figs. 16b and 16d, for rriA — ► there is a positive lower bound on sm 2 (/3 — a) 
(see fig. 7), reflecting the fact that at the tree level a — > —(3 in this limit, 
so the SM value is actually never reached. For increasing m#, one can 
notice the fast decoupling of H from W-pairs, as already observed when 
discussing the total width. In leading order, A does not couple to iy-pairs. 
A non- vanishing cross-section could be generated at one loop, but such a 
contribution is completely negligible, since even for h and H the W-fusion 
cross-section is small (^ 20%) compared with the gluon-fusion cross-section. 
Finally, we observe that the LHC and SSC cross-sections of fig. 16 differ by 
an overall factor ~ 3 in the phenomenologically relevant region, = 70-140 
GeV. 
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Associated bb(p production. This mechanism is unimportant in the 
SM, since its cross-section is too small to give detectable signals [Q. In 
the MSSM model, however, for large values of tan/3 the bbcj) couplings can 
be strongly enhanced. Then for not too high values of the Higgs masses, a 
significant fraction of the total cross-section for neutral Higgs bosons can be 
due to this mechanism. 

The associated bb<j> production involves two rather different mass scales, 
m<£ >> rut, therefore at higher orders large logarithmic corrections of order 

a>" 

\mhJ 

may destroy the validity of the Born approximation, depending on the value 
of m<f). One needs an improved treatment where these logarithms are re- 
summed to all orders. The origin of these logarithms is well understood. Part 
of them come from configurations where the gluons are radiated collinearly 
by nearly on-shell bottom quarks, which are obtained by splitting the initial 
gluons into a bb pair. This type of logarithms are responsible for the QCD 
evolution of the effective 6-quark density within the proton: they were care- 
fully analysed and resummed to all orders, and it was found that the 
corrections are positive and increase with the Higgs-boson mass. A second 
subset of logarithms lead to running quark mass effects. An analysis where 
both effects are treated simultaneously is still missing. In view of this ambi- 
guity, we interpolated the existing results by using the Born approximation 
with the bottom quark mass adjusted to the fixed value m& = 4 GeV. How- 
ever, one should keep in mind that the theoretical estimate in this case has 
a large (factor of 2) uncertainty. 

In fig. 17 we display cross-sections for associated bbcj) production, for the 
same tan/3 and y/s values as in fig. 15. Comparing the cross-sections of 
figs. 15 and 17, we can see that the hbb cross-section can give at most a 20% 
correction to the total h cross-section. The Hbb and Abb cross-sections, how- 
ever, can be even larger than the corresponding gluon-fusion cross-sections 
for tan (3 10. Comparing the LHC and SSC curves of fig. 17, one can notice 
a rescaling factor varying from ~ 3 to ~ 8 in the region from 60 to 500 
GeV. 

Associated W<p (Z<j>) production. This mechanism 0] is the hadron 
collider analogue of the SM Higgs production mechanism at LEP, with the 
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difference that at hadron colliders Wcp production is more important than 
Zip production. In the Zip case, the event rate at the LHC and SSC is 
too low to give a detectable signal, both in the SM and (consequently) in 
the MSSM. The W<p mechanism has considerable importance at the LHC 
for <p = h,H and in the Higgs mass range = 70-140 GeV, where a 
measurable signal may be obtained from final states consisting of two isolated 
photons and one isolated lepton. The calculation of the cross-section is well 
understood, including the QCD corrections, since it has a structure similar to 
the Drell-Yan process, with the same next-to-leading-order corrections (for 
a recent study concerning the numerical importance of the QCD corrections 
see ref. f4l[]). The QCD corrections are positive, and amount to about 12% 
if one chooses Q 2 = s as the scale of Q 2 evolution. The production cross- 
sections of h and H are obtained by rescaling the SM model cross-section by 
the appropriate correction factors given in table 1. 

In fig. 18, cross-sections for Wh and WH are displayed, as functions of 
corresponding Higgs masses, for the same tan/3 and y/s values as in fig. 15. 
Since the SUSY correction factors are the same, the approach to the SM case 
and the irrelevance of WA production can be described in the same way as 
for the W^-fusion mechanism. 

In the phenomenologically relevant region, = 70-140 GeV, the scaling 
factor between the LHC and SSC curves is ~ 2.5. 



Associated ti<p production. In the SM, the Born cross-section 
formula for this process is the same as for the bbip case |[38| . In the MSSM 
case, one just needs to insert the appropriate SUSY correction factors, as 
from table 1. Note, however, that the leading-order QCD calculation is more 
reliable in this case, since in the ticj) case one does not have two very different 
physical scales when is in the intermediate mass region. The next-to- 
leading QCD corrections are not known, therefore the Born cross-section 
still suffers from a relatively large (~ 50%) scale ambiguity. 

In fig. 19, the production cross-sections for ticj) (</> = h, H, A) are plotted, 
as functions of the corresponding Higgs mass, for the same tan/5 and a/s 
values as in fig. 15. In general, the MSSM cross-sections are smaller than the 
SM one, which is approached in the limit in which the ti<p coupling becomes 
SM-like. A possible exception is the ttH cross-section for small values of 
rriA and tan/?, since in this case the corresponding coupling can be slightly 
enhanced with respect to the SM one. 
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In the phenomenologically relevant range, = 70-140 GeV, the rescal- 
ing factor between the LHC and SSC curves in fig. 19 varies from ~ 6 to 
~ 7. 



In the phenomenologically allowed range of eq. ([1]), the top- mass depen- 
dence of the cross-sections of figs. 15-19 is not negligible, but it does not 
change qualitatively the previous considerations. The largest effect comes 
from the increase of the upper limit on rrih for increasing top mass (see 
fig. 1). This induces a shift in the limiting values for the h and H production 
cross-sections. There are also obvious kinematical top-mass effects in the 
gluon-fusion mechanism and in the titfi mechanism, which are well under- 



stood from SM studies Ul4| . In the MSSM, additional effects are given by the 
radiative corrections to the relevant Higgs couplings, which were discussed 
in sect. 2. 



6 Physics signals 

6.1 Neutral Higgs bosons 

We now calculate the rates for a number of processes that could provide 
evidence for one or more of the neutral MSSM Higgs bosons at the LHC 
and SSC, and we summarize our results with the help of contour plots in 
the (mA,tan/3) plane. We consider production cross-sections, folded with 
branching ratios, for the following signals: 

• two isolated photons; 

• one isolated lepton and two isolated photons; 

• four isolated charged leptons; 

• a pair of tau leptons. 

In the SM case, the first two signals are relevant for the region of intermediate 
Higgs mass, 70 GeV < < 140 GeV, the third one is the so-called 'gold- 
plated' signal in the high-mass region 130 GeV < < 800 GeV, and the 
t + t - signal appears to be hopelessly difficult. 
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In a complete phenomenological study, one would like to determine pre- 
cisely the statistical significance of the different physics signals. This would 
require, besides the computation of total signal rates, the calculation of the 
backgrounds, the determination of the efficiencies (for both signals and back- 
grounds) due to kinematical cuts and detector effects, the optimization of 
the kinematical cuts to achieve the best signal/background ratio, etc. Such 
a complete analysis would require the specification of several detector and 
machine parameters, and goes beyond the aim of the present paper. Instead, 
we try here to present total rates for well-defined physics signals, in a form 
which should be useful as a starting point for dedicated experimental studies. 

As the only exception, to illustrate with an example how our results can 
be used to establish the statistical significance of a given physics signal in a 
given detector, we shall describe the case of the 'two-isolated-photons' signal, 
using the results of recent simulation works. A similar procedure should be 
adopted for any other physics signal, detector, and collider, once complete 
results of simulation works are available. In many cases, the existing results 
from previous background and simulation studies, carried out for the SM, can 
also be used to draw conclusions concerning the MSSM case. We mention, 
however, two important differences: 1) the total widths of H and A remain 
small even in the high-mass region; 2) for large tan/3, the number of signal 
events in the r + r~ final state is significantly higher than in the SM case. 



Inclusive two-photon channel. In fig. 20 we display cross-sections 
times branching ratios for the inclusive production of <fi — h, H, A, followed 
by the decay <p —> 77; as functions of m^, and for the same parameter choices 
and energies as in fig. 15. We sum the contributions of the gluon-fusion, In- 
fusion, and bb(f) mechanisms. For comparison, the SM value is also indicated. 



The QCD corrections of ref. |3TJ] are not included, for the reasons explained 



in the previous section. In the case of h and H, the signal rates are always 
smaller than in the SM, and approach the SM values at the upper and lower 
edge of the allowed m h amd m H ranges, respectively. The rather steep slope 
characterizing the approach to the SM limit, for varying Higgs mass and 
fixed tan /3, is a reflection of the property of the branching ratios discussed 
in sect. 4. Also the structure in figs. 20b and 20e can be attributed to the 
threshold behaviour of the H —>■ hh channel. The signal rate for the CP-odd 
A boson is extremely small for tan/5 <; 3. However, we observe that in a 
small region of the parameter space, for rriA just below 2m t and tan/? ^ 3, 
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the rate can become larger than the SM one: nevertheless, in general it is 
still too low to produce a detectable signal, unless one chooses tan (3 ~ 1 and 
rriA ~ 2m t . 

For the inclusive two-photon channel, the results of detailed simulations 
of signal and SM background are now available, for some of the LHC detector 



concepts [42,16,43,44]. For illustrative purposes, in the following example we 



shall follow the treatment of ref. fl43l . In the mass range = 80-150 GeV, 
and assuming 10 5 pb _1 integrated luminosity, this LHC simulation considers 
a fairly wide range of detector performances, which affect the significance 



of the signal. For an energy resolution AE/E = [2%/^ E(GeV)] + 0.5%, 
ref. [fKJ obtains a 10 4 efficiency for rejecting jets faking an isolated photon 
in the relevant px region. Applying standard kinematical cuts, this simula- 
tion finds ~ 40-50% kinematical acceptance, with an additional ~ 30-40% 
loss due to isolation cuts and reconstruction efficiency for the isolated pho- 
tons. Typically, for a SM Higgs with ~ 100 GeV, one obtains ~ 10 3 
signal events over ~ 10 4 background events, corresponding to a statistical 
significance S/\fB ~ 10. More generally, ref. |^3| determined the statistical 
significance of the signal for given values of the generic Higgs mass and 
of the signal rate a ■ BR(<fi — > 77) (see fig. 21). In our opinion, this is an ex- 
cellent way of summarizing the simulation work, since it gives the possibility 
of studying alternatives to the SM case, and in particular the MSSM. The 
dashed line in fig. 21 corresponds to the signal for the SM Higgs, which in- 
cludes both the gluon-fusion and the H^-fusion production mechanisms, and 



also the QCD corrections of ref. [ 31| . One can see from fig. 21 that for such 
optimistic detector parameters there is some margin for detecting smaller 
rates than in the SM. Clearly the SUSY Higgs search further enhances the 
need for the best possible m 77 resolution and 7-jet rejection. 

In extending the SM analysis to the MSSM, one should pay attention 
to the applicability of the QCD corrections of ref. |3l| to the gluon-fusion 
cross-section. We have checked that in the phenomenologically relevant re- 
gion, which corresponds to h or H in the intermediate mass range, and to 
signal rates within an order of magnitude from the SM one, the gluon-fusion 
mechanism is dominated by the top-quark loop. Since in this region the 
gluon-fusion mechanism accounts for ~ 80% of the total cross-section, and 
the correction is roughly a multiplicative factor 1.5, as a rule of thumb we 
can take it into account by multiplying the total cross-section by a factor 
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~ 1.4. 

In fig. 22 we show contour plots in the (tua, tan f3) plane, corresponding 
to fixed values of a ■ BR((p — > 77) (</> = h, H). QCD corrections have been 
included according to ref. [[H]]. The region where the rate is large enough to 
promise a measurable signal is rather large for h, is concentrated in a small 
strip for H, and is possibly a very small area, just below = 2m t and 
just above tan (3 = 1, for A. In our representative example [43], we can now 
evaluate the statistical significance of the 'two-isolated-photons' signal at any 
point of the (m^,tan/3) plane, by just combining the information contained 
in figs. 22, 2, 3, and 21. In the case of h searches, and for ^ 90 GeV, a 
signal rate beyond 40 fb should give detectable signals. A signal rate of 30 
fb is the borderline of detectability for one year of running, and signal rates 
below 20 fb appear extremely difficult to detect. In the case of H, which has 
higher mass, a signal rate of 20 fb appears to be the borderline of what can 
be achieved in one year of running. In the case of A, the interesting mass 
region is itla ~ 2m t : for rriA = 250 GeV, and taking a ■ BR(A 77) = 3 fb 
as a plausible discovery limit at the LHC f4"3]| , a signal for A — ► 77 will be 
found only if tan/5 ^ 1.5. 

One isolated lepton and two isolated photons. This signal can 
come from either Wcj) or ti(f) production. In the latter case, two or more 
isolated jets are also produced. The physics signals from Wcp production are 
particularly important at the LHC, and were studied in ref. fl4"5"f . The im- 



portance of the physics signals from ttcp production was recently emphasized 



in ref. ||46|| . The production rates, multiplied by the <p — > 77 branching ratio, 
are shown in fig. 23. We can see that, similarly to the inclusive 77 channel, 
the rates for Wh, WH, tih, tiH are always smaller than the SM value, which 
represents the boundary curve in the limit of large tan/3 and large (small) 
rriA for h (H). From figs. 23g and 23j one can see that it A production can 
give a /Z7 signal larger than in the SM for small tan (3 and near the tjia = 2m t 
threshold, but even in this case the rate appears to be too small for detection. 
We emphasize that the production rates shown in fig. 23 do not include the 
branching ratios of leptonic W and semileptonic t decays. If top decays are 
as in the SM, one should still include a combinatorial factor of 2, coming 
from the fact that both top and antitop can decay semileptonically. On the 
other hand, in the MSSM there is the possibility of t — > bH + decays, where 
the subsequent H + decay cannot produce a direct lepton I — e, We shall 
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take this possibility into account in the following, but its impact on the de- 
tect ability of the £77 signal is rather small. The only case in which this effect 
is not completely negligible is for H, when ^ 100 GeV and tan (5 ^ 4 or 
tan f3 ^ 10, in which case the t — > bH + branching ratio can play a role. 

From parton- level simulations |fi2l , [47 t fl8| , for a SM Higgs of about 100 
GeV, one typically obtains ~ (12 + 15) and ~ (3 + 11) Z77 signal events 
at the LHC and SSC, respectively. Here we assumed 10 5 pb _1 of integrated 
luminosity for the LHC and 10 4 pb _1 for the SSC. The quoted numbers sep- 
arately show the contributions from Wp> and tiip production. Furthermore, 
they include losses due to acceptances (~ 30%), and lepton and photon de- 
tection efficiencies [e ~ (0.9) 3 ]. The total background is roughly 20-30% 
of the signal and is dominated by the irreducible ^77 and ££77 contribu- 
tions. There are many different contributions to the reducible background 
(bbg, bby, 6677, Wjj, . . .). Parton-level simulations indicate that they can be 
suppressed well below the irreducible background, provided one assumes, as 
for the inclusive 77 case, excellent detector performances: a 7-jet rejection 
factor ^ 3 x 10 3 and a suppression factor <: 7 for the leptons from 6-decays 
after isolation cuts. 

Clearly, there is very little margin (a factor of 2?) to be sensitive to signal 
rates smaller than in the SM. In fig. 24, we show contour plots corresponding 
to fixed values for the quantity 

= [2a{U4>) ■ BR(t -> Wb) + a(W<j>)] ■ BR(<j) -> 77) ■ BR(W -> Iv) , 

(40) 

for the same choice of m t and rriq as in fig. 15, and for LHC and SSC energies. 
In eq. (f40"D, I = e, /x and we have not considered the strongly suppressed 
possibility of getting a light charged lepton from both top and antitop. 

The four-lepton channel. The channel ip — > Z*Z* — > /+/ /+/ 
(7 = e, fx) gives the so-called 'gold-plated' signal for the SM Higgs in the 
mass range = 130-800 GeV. Below m v ~ 130 GeV, both the total rate 
and the acceptance decrease very rapidly, leading to too small a signal for 
detection. For all three neutral Higgs bosons of the MSSM, the rates in this 
channnel are always smaller than in the SM. In the case of A, there is no 
AZZ coupling at tree level, and loop corrections cannot generate measurable 
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rates in the four-lepton channel. As for h, if m t ^ 180 GeV and rriq ^ 1 TeV, 
one can see from fig. 1 that m^ ^ 130 GeV. Therefore, the h — > Z*Z* — > 41 
signal does not have chances of detection at the LHC and SSC, unless one 
chooses extremely high values for m t and rriq or one has superb resolution 
and acceptance for leptons. The situation is somewhat better in the case 
of H, despite the strong suppression with respect to the SM, due to the 
competition with the hh,bb,tt channels, as discussed in sect. 3. 

In fig. 25, we show signal rates for the SM Higgs boson and for H, for 
the same choices of parameters as in fig. 15. The threshold effects and the 
suppression for large values of tan /3 are clearly visible. 

The LHC and SSC discovery potential can be estimated by using the 
results of simulations carried out for the SM P^|r45|j5Tl|JTB|JTE[ ] , taking also into 



account that Th < 2 GeV all over the mass region of interest, m# ^ 2m t . 
Assuming excellent lepton momentum resolution, in the mass range 2m>z < 
mu < 2m t a signal rate ~ 20 smaller than in the SM could still lead to a 
detectable signal. In fig. 26, we show contour plots in the (m^, tan (3) plane, 
corresponding to fixed values of a ■ BR(H — > 4/). QCD corrections have been 



included according to ref . |3T| . In view of the strong sensitivity to the value of 
m t , we show contours for m t = 120, 140, 160 GeV, for LHC and SSC energies, 
and for rriq = 1 TeV. The two almost vertical dashed lines correspond to 
mu = 2m z and to m# = 2m t . For mu > 2m z, a detectable signal could be 
obtained up to tan/3 ~ 5. Notice that the experimental acceptances change 
with mn\ in particular, in the region mu < 2m z they fastly decrease with 
decreasing mu'- for a realistic assessment of the discovery limits in this mass 
region, one should take this and other effects into account. Anyway, the 
prospects for detection for mu < 2mz do not look good if m t ^ 150 GeV 
and rriq ^ 1 TeV. 

The t + t~ channel. For the SM Higgs boson, the r + r~ decay channel 
has been found hopelessly difficult for discovery |51 , 52] , since this channel has 



bad mass resolution and overwhelmingly large background coming from the 
production of ti, WW+ jets, Drell-Yan pairs, Z+ jets, bb+ jets, .... The bad 
resolution is due to the fact that the tau-decay products always include one or 
more neutrinos, which carry away energy; therefore one cannot reconstruct 
the signal as a resonance peak. The situation is improved if the Higgs is 



produced with large transverse momentum that is balanced by a jet [53 . In 
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this case one can use the approximation 



-(1) -(2) 
(1) , (2) Pl£_ _ -jmiM / 41 x 

(i) + ( 2 ) — Pi K^ 1 -) 

to reconstruct the transverse momenta of the neutrinos and hence the in- 
variant mass of the tau pair. In the above equation, pW denotes the total 
transverse momentum of the neutrinos coming from the decay of i — 1, 2, 
while Pij, and £j denote the lepton momenta and energies, respectively. It 
was shown in ref. that, in the mass range = 70-140 GeV, a mass 
resolution of ~13— 17% can be achieved. This method can also be used for the 
hadronic decay modes, taking advantage of the fact that the rate is higher 
by a factor of ~ 3.5. When a tau decays hadronically, the hadrons have very 
low multiplicity and invariant mass, and these properties might be used to 



recognize the V-jet' There is a price for the better mass resolution. 

Tagging on a large-px jet can reduce the rate by an order of magnitude. Fur- 
thermore, at 10 5 pb~ /year luminosity, the presence of pile-up deteriorates 
significantly the measurement of Pt^, and therefore the r + r~ signal can 
only be studied with this method at lower luminosities, ~ 10 4 pb _1 /year Q. 
While these difficulties appear prohibitive in the case of the SM, the situation 
is not entirely negative in the MSSM. 

In fig. 27, we display signal rates for </> — > t + t~ production (0 = h, H, A), 
for the same parameter choices as in fig. 15, together with the SM values. 
We can see that for large values of tan/3 the production rates can become 
much larger than in the SM. In the case of h production, for tan/3 = 30 
the enhancement can be more than one order of magnitude, and increases 
with decreasing values of rrih (see figs. 27a and 27d). Huge enhancements 
can be obtained also for H and A, thanks to the properties of the t + t~ 
branching ratios discussed in sect. 4. Note for example that, at the LHC, for 
m H , m A ~ 500 GeV and tan/3 > 10, we get a ■ BR(H, A -> r+r - ) > 20 pb, 
while for rriA^H ~ 120 GeV and tan/3 ^ 30, we obtain a ■ BR(H,A — ► 
t + t - ) > 30 pb. The rates for the SSC are rescaled by the factor already 
discussed in the previous section. 

In order to assess, for any given mass, the cross-section values above 
which one obtains a measurable signal over the large background, detailed 

Alternatively, at high luminosity one may try to just search for an excess of events in 
the e jj^ or V-jet' channels. 
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simulations are needed. Preliminary studies have been reported for the lepton 
channel e ± /i T in ref. [|5l] and for the mixed channel / ± + 'r-jet' in ref. 



In the second case, the difficulty of recognizing a 'r-jet' may be compensated 
by the higher rate of this channel. The preliminary analysis of ref. |55j finds 



for the LHC sensitivity to values of a ■ BR((p — > t + t~) down to ~ 10 pb 
in the low-mass region ~ 100 GeV and ~ 1 pb in the high-mass region 
m<£ ~ 400 GeV. This result cannot be easily rescaled to the SSC case, since a 
large mass interval is involved and the SSC luminosity gives more favourable 
experimental conditions for the srudy of this channel. 

In fig. 28, we show contour plots corresponding to fixed values of o ■ 
BR(<p — > t + t~), for the same values of m t and rriq as in fig. 15. 

6.2 Charged Higgs boson 

We now move to the discussion of possible physics signals associated with 
the charged Higgs boson. The phenomenology of the charged Higgs boson 
at hadron colliders was previously discussed in refs. The benchmark 



mass value for charged-Higgs-boson searches at the LHC and SSC is m#± = 
m t — mj. For lower values of ttih±, the dominant production mechanism 
at large hadron colliders is gg — > it, followed by t — > H + b. For higher 
values of ra#±, the dominant production mechanism is gb — > tH + . As far as 
detectable signals are concerned, this last case appears hopeless, in view of 
the suppressed cross-section and of the large backgrounds coming from QCD 
subprocesses. The first case appears instead to be experimentally viable over 
a non-negligible region of parameter space. Given the known it production 
cross-section, one can compute the t — > bH + branching ratio according to 
well-known formulae, parametrizing again the leading QCD corrections by 
running masses evaluated at a scale Q ~ m H ±. The charged Higgs branching 
ratios were discussed in the previous section, where it was found that the t + u t 
mode dominates in the mass range under consideration. The experimental 
signal of a charged Higgs would then be a violation of lepton universality 
in semileptonic top decays. As a convenient indicator, one can consider the 
ratio 

Bfi (^T + „ r fr) 

BR(t -> n+v^b) ' V ' 
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with 

BR{t -> H+b) ■ BR(H + -» r+z/ T ) 

fli2(t -> • -> //+!/„) ' 1 ' 

Preliminary investigations fHfj] show that the experimental sensitivity could 
reach Ai? ^ 0.15 at the LHC. At the SSC the increased tt production cross- 
section is likely to give better sensitivity In fig. 29, we display contour 
lines of AR in the (rriA, tan 0) plane, for the three representative values 
m t = 120, 140, 160 GeV. The dashed lines denote the kinematical limit 
m,H± = m t — mj. One can see that the most difficult values of tan (3 are those 
between 3 and 10, and that the process under consideration could give access 
to values of itia as high as 80-120 GeV for top-quark masses in the range 
120-160 GeV. 



7 Conclusions and outlook 

In this paper we carried out a systematic analysis of the possible physics 
signals of the MSSM Higgs sector at the LHC and SSC, assuming that the 
super symmetric (R-odd) particles are heavy enough not to affect significantly 
the production cross-sections and the branching ratios of the MSSM Higgs 
particles. As independent parameters in the Higgs sector, we chose rriA and 
tan p, and we considered the theoretically motivated region of the parameter 
space 

< m A < 500 GeV, 1 < tan/3 < 50 . 

We assumed rriq = 1 TeV and negligible mixing in the squark sector. We 
included the most important radiative corrections to the Higgs masses rrih, 
m H, Wiff±, and to the Higgs couplings to fermions and vector bosons. We 
also included the most important radiative corrections to the three-point 
couplings of the neutral Higgs bosons. 

We estimated the discovery potential of LEP I and LEP II, and we carried 
out detailed cross-section calculations for the LHC and SSC. We singled out 
four classes of final states (77, / ± 77, r + r _ ) which could provide 

significant signals for neutral Higgs bosons at the LHC and SSC, and we also 
examined possible signals of charged Higgs bosons in top decays. 

We calculated all the relevant branching ratios, and the cross-sections for 
all the relevant production mechanisms. We presented our results with the 
help of branching-ratio curves, cross-section curves, signal-rate curves and 
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contour plots in the (mi, tan (3) plane. We did not perform new background 
studies, but we pointed out that, using the results of our calculations and 
of the existing simulations carried out for the SM Higgs, supplemented by 
estimates of the acceptances and efficiencies of typical experiments, in many 
cases one can draw conclusions concerning the discovery range. In some cases, 
as for the r + r~ channel, further simulation work appears to be needed in 
order to reach firm conclusions. Nevertheless, some preliminary conclusions 
can already be drawn and will now be summarized. 

At large hadron colliders, the MSSM Higgs search is, in general, more 
difficult than the SM Higgs search. This is due to the fact that, in a large 
region of the parameter space, at least one of the MSSM neutral Higgs bosons 
is in the intermediate mass region, 80 GeV ^ 140 GeV, but with rates 

in the 77 and Z77 channels which can be significantly suppressed with respect 
to the SM case. Similarly, neutral Higgs bosons with ^ 130 GeV have 
typically strongly suppressed rates in the channel. On the contrary, 

in the MSSM, for rather large values of tan /3, one can obtain a much larger 
signal rate in the r + r~ channel than in the SM. Finally, t — > bH + decays, 
followed by H + — > t + u t , can give detectable signals only in a rather restricted 
region of the parameter space. 

As an example, we now try to assess the discovery potential of the different 
channels for the representative parameter choice m t = 140 GeV, trig = 1 TeV, 
working as usual in the (m^tan/?) plane. 

To begin with, we recall the expectations for LEP II. The size of the 
LEP II discovery region depends rather strongly on m t and m,g, and on the 
assumed energy and luminosity. For standard machine parameters, LEP II 
cannot test the whole parameter space allowed by the present data. Looking 
back at fig. 8, one may tentatively say that LEP II will give us (if no Higgs 
boson is discovered) lower limits of about tua *J 70-100 GeV and tan/3 
3-8 for m t = 120 GeV, tan/3 £ 1.5-3 for m t = 160 GeV. 

We observe that the LHC and SSC will test the MSSM Higgs sector in a 
largely complementary region of the (rattan/?) plane. A pictorial summary 
of the discovery potential of the different channels is presented in fig. 30. 
We emphasize once again that the final discovery limits will depend on the 
machine and detector properties, as well as on the actual values of the top 
and the soft supersymmetry-breaking masses. We therefore drew fig. 30 
just for illustrative purposes, to exemplify a particularly convenient way of 
considering all the discovery channels at once. 
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The 77 and Z77 channels are important in approximately the same region 
of the parameter space, ttia 3> rriz for h and 50 GeV ^ uia ^100 GeV 
for H. Therefore, these two channels can be experimentally cross-checked 
one against the other, reinforcing the significance of a possible signal. As an 
optimistic discovery limit for h, we show in fig. 30 the contour line <j-BR(h — > 
77) = 30 fb at the LHC, corresponding to rriA <; 200 GeV. This contour 
line is shown only for irih ^ 80 GeV. In the region of the parameter space 
to the right of this line [indicated by the labels h — > 77 and / + (/&—> 77)], 
it is expected that measurable signals will be found, assuming detector and 
machine parameters as discussed in refs. |16|]. Approximately the same 
contour line is obtained by taking a ■ BR{h — > 77) = 85 fb at the SSC. This 
indicates that, in the inclusive 77 channel, the discovery range of the LHC 
and SSC will be the same if the luminosity at the LHC will be ~ 3 times 
higher than at the SSC and if the detectors used at the two machines will have 
similar efficiencies in suppressing the backgrounds. Very similar discovery 
lines can be drawn by considering the Z77 channel and taking o ■ BR[l + {h — > 
77)] ~ 0.8 fb for the LHC and ~ 4 fb for the SSC. 

In fig. 30 we also show the contour line for a ■ BR(H — > 77) = 20 fb at 
the LHC, corresponding to ~ 55 fb at the SSC. The slightly smaller signal 
rate was chosen to account for the improved efficiencies at higher Higgs-mass 
values. The contour line defines a narrow strip around rriA ~ 75 GeV (shaded 
in fig. 30), where the discovery of H is expected to be possible both in the 
77 and in the £77 channels [for lack of space the label l + (H —> 77) has been 
omitted] . 

The four-lepton channel is important mainly for H, in the mass region 
2m z ^ ran ~ 2m t , which translates into 150 GeV ^ rriA ~ 2m t , and 
for relatively small tan (3. As a reference value for discovery in this mass 
region, we take o ■ BR(H —>■ 4/) = 1 fb for the LHC, which corresponds to 
a ■ BR(H — > 4/) ~ 3 fb for the SSC. This contour defines the area in fig. 30 
indicated by the label H — > 41. In a small part of this area, corresponding 
to niA ~ rat and tan/3 ~ 1, A — > 77 could also give a detectable signal. 

In the region of very large tan f3, and moderately large tua, one could take 
advantage of the enhanced production cross-sections and of the unsuppressed 
decays into r + r _ to obtain a visible signal for one or more of the MSSM 
neutral Higgs bosons, and in particular for H and A, whose masses can be 
significantly larger than 100 GeV. The simulation work for this process is 



still at a rather early stage [55], so that no definite conclusion can be drawn 
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yet. For reference, the dotted line in fig. 30 corresponds to a (somewhat 
arbitrary) interpolation of a ■ BR(<f) — > t + t~) ~ 10 pb at = 100 GeV and 
a ■ BR((j) — > r + r~) ~ 1 pb at = 400 GeV, for LHC energy and summing 
over the </> = H,A channels. 

Finally, in the region of parameter space corresponding to ^100 GeV, 
a violation of lepton universality due to the decay chain t — ► bH + — > br + u T 
could indicate the existence of the charged Higgs boson. This region is in- 
dicated by the label H^ — > rv in fig. 30. Its right border is defined by the 
contour line of R = 1.15, where R was defined in eq. (f42"l). 

By definition, our contour lines do not take into account changes in the 
acceptances and efficiencies, which are expected in realistic experimental con- 
ditions, and depend on the Higgs mass and on the detector. We therefore 
expect some deformations of our contours once discovery lines are extracted 
from realistic experimental simulations |33],55]]. 



As a last piece of information we also display in fig. 30 the border of 
the expected discovery region at LEP II, which depends rather sensitively, as 
already discussed, on the assumed values of the machine energy and luminos- 
ity. We then show two representative lines: the lower dashed line corresponds 
to a(e + e- -> hZ, HZ, hA, HA) = 0.2 pb at ^ = 175 GeV, while the up- 
per dashed line corresponds to cr(e + e _ — > hZ, HZ, hA, HA) = 0.05 pb at 
y/s = 190 GeV. 

Using the result summarized in fig. 30, we can draw several important 
qualitative conclusions: 

• The discovery potentials of LEP and the LHC/SSC show a certain 
complementarity. The discovery region at LEP covers all tan (3 values 
at small values of uia, and all uia values at small values of tan f3, while 
at the LHC/SSC one should be sensitive to the large tan/3, large 
region. 

• One may ask whether the LHC and SSC, combined with LEP II, can 
explore the full parameter space of the MSSM Higgs sector, being sen- 
sitive to at least one signal in each point of the (m^, tan/3) plane, for 
all plausible values of mt and rriq. At present, this question cannot be 
answered positively. The union of the regions where one should find 
signals at least for one Higgs boson does not cover the whole parameter 
space: the discovery region has a hole in the middle of the parameter 
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space. For our parameter choice, the most difficult region appears to 
be the cross-hatched area around uia = 150 GeV and tan/5 = 10. 
Therefore we cannot claim yet the existence of a 'no- lose' theorem for 
the MSSM Higgs search. 

• One may also ask if there are regions of parameter space where one can 
find more than one signal from the MSSM Higgs sector. The answer 
is that around = 200 GeV and tan (3 < 5 one can discover h at 
LEP II and H at the LHC/SSC in the four-lepton channel. There is 
a somewhat smaller region above m\ = 200 GeV where one can also 
find h in the 77 and Z77 channels. Furthermore, at high values of tan (3 
(<: 20) and at tua > 200 GeV one may discover A and H in the tt 
channel and h in the 77 and Z77 channels, although the separation of H 
and A appears to be impossible, due to their almost perfect degeneracy 
in mass. This part of the parameter space is inaccessible at LEP II. The 
discovery region for H in the 77 and £,77 channels, corresponding to 
low values of m A , largely overlaps with the LEP II discovery region and 
with the discovery region related to charged-Higgs production in top 
decays. In the low tan/3, for 80 GeV ^ ^ 180 GeV and ^ 2mi 
one has a signal at LEP II and no signal at the LHC and SSC, since 
rrih is too small for detection. 

Finally, we would like to make some comments on the theoretical uncer- 
tainties and on possible future studies. 

Our values for the signal rates depend on several phenomenological input 
parameters, as the value of the bottom mass, the parton-number densities 
and the value of as- The given cross-sections and branching ratios will 
change if the input parameters are varied in their allowed range. Also, for 
some production mechanisms, only the Born cross-sections are known. We 
estimate that the theoretical errors of the calculated rates vary from about 
30%, in the case of the 77 channel, up to about a factor of 2 when the bbcf) 
or tt(j) production mechanisms are important. 

We did not study all effects associated to variations of the parameters in 
the SUSY (i?-odd) sector. It would be interesting to examine the case when 
some of the Higgs bosons are allowed to decay into i?-odd SUSY particles, 
or the effects of squark mixing. More importantly, serious simulation work 
is still needed, in particular for the r + r~ and the Z77 channel. 
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Note added 



After the completion of most of the work presented in this paper, which was 



anticipated in many talks ||57|| , we received a number of papers [|8|,|)£| dealing 
with different subsets of the material presented here, and reaching similar 
conclusions. Reference |5S] also contains the generalization of eqs. (|21~| -f27|) 
to the case of arbitrary mixing in the stop and sbottom mass matrices, but 
still neglecting the D-term contributions. 
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Figure captions 



Fig.l: Contours of m™ ax (the maximum value of m h , reached for m A — > oo): 

a) in the (m t ,tan/3) plane, for rriq = 1 TeV; b) in the (m t ,mq) plane, 
for tan/3 = m t /mb. 

Fig.2: Contours of m h in the (m A ,ta.n (3) plane, for rriq = 1 TeV and a) m t = 
120 GeV, b) m t = 160 GeV. 

Fig.3: Contours of ran in the (rattan/?) plane, for rriq = 1 TeV and a) m t = 
120 GeV, b) m t = 160 GeV. 

Fig.4: Contours of m H ± in the (m A ,t&n(3) plane, for rriq = 1 TeV. The solid 
lines correspond to m t = 120 GeV, the dashed ones to m t = 160 GeV. 

Fig.5: Contours of sin 2 a/ cos 2 /3 in the (m^tan (3) plane, for rriq = 1 TeV. 
The solid lines correspond to m t = 120 GeV, the dashed ones to m t = 
160 GeV. 

Fig.6: Contours of cos 2 a/ sin 2 /3 in the (m,4,tan (3) plane, for rriq = 1 TeV. 
The solid lines correspond to m t = 120 GeV, the dashed ones to m t = 
160 GeV. 

Fig. 7: Contours of sin 2 (/9 — a) in the (m A ,tan(3) plane, for rriq = 1 TeV. 
The solid lines correspond to m t = 120 GeV, the dashed ones to m t = 
160 GeV. 

Fig.8: Schematic representation of the present LEP I limits and of the fu- 
ture LEP II sensitivity in the (m^tan (3) plane, for rriq = 1 TeV 
and a) m t = 120 GeV, b) m t = 160 GeV. The solid lines corre- 
spond to the present LEP I limits. The dashed lines correspond to 
a(e + e- -> hZ, HZ, hA, HA) = 0.2 pb at ^1 = 175 GeV, which could 
be seen as a rather conservative estimate of the LEP II sensitivity. The 
dash-dotted lines correspond to cr(e + e~ — > hZ, HZ, hA, HA) = 0.05 pb 
at t/s = 190 GeV, which could be seen as a rather optimistic estimate 
of the LEP II sensitivity. 

Fig.9: Total widths of the MSSM Higgs bosons, as functions of their respective 
masses, for m t = 140 GeV, rriq = 1 TeV and tan (3 — 1.5, 3, 10, 30: a) h; 

b) H; c) A; d) H ± . 
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Fig. 10: Branching ratios for h, as functions of nth, for m t = 140 GeV, rriq = 
1 TeV and: a) tan/5 = 1.5; b) tan/3 = 30. 

Fig. 11: Branching ratios for H, as functions of m#, for m t = 140 GeV, rriq = 
1 TeV and: a) tan (5 = 1.5; b) tan (5 = 30. 

Fig. 12: Branching ratios for A, as functions of m^, for m t = 140 GeV, rriq = 
1 TeV and: a) tan (5 = 1.5; b) tan/3 = 30. 

Fig. 13: Branching ratios for h, as a function of tua, for m t = 140 GeV, rriq = 
1 TeV and: a) tan (5 = 1.5; b) tan (5 = 30. 

Fig. 14: Branching ratios for H ± , as functions of m H ±, for m t = 140 GeV, 
rriq = 1 TeV and: a) tan/3 = 1.5; b) tan/3 = 30. 

Fig. 15: Cross-sections for neutral Higgs production, via the gluon-fusion mech- 
anism, as functions of the corresponding masses and for m t = 140 GeV, 
rriq = 1 TeV, tan/3 = 1.5,3, 10,30: a) h, LHC; b) H, LHC; c) A, LHC; 
d) h, SSC; e) H, SSC; f) A, SSC. QCD corrections are not included. 

Fig. 16: Cross-sections for h and H production, via the VT-fusion mechanism, as 
functions of the corresponding masses, for the same parameter choices 
as in fig. 15: a) h, LHC; b) H, LHC; c) h, SSC; d) H, SSC. 

Fig. 17: Cross-sections for associated btxj) production, as functions of the corre- 
sponding Higgs masses, for the same parameter choices as in fig. 15: 
a) h, LHC; b) H, LHC; c) A, LHC; d) h, SSC; e) H, SSC; f) A, SSC. 

Fig. 18: Cross-sections for associated W(j) production, as functions of the corre- 
sponding Higgs masses and for the same parameter choices as in fig. 15: 
a) h, LHC; b) H, LHC; c) h, SSC; d) H, SSC. 

Fig. 19: Cross-sections for associated ticj) production, as functions of the corre- 
sponding Higgs masses and for the same parameter choices as in fig. 15: 
a) h, LHC; b) H, LHC; c) A, LHC; d) h, SSC; e) H, SSC; f) A, SSC. 

Fig. 20: Cross-sections times branching ratios for inclusive Higgs production 
(the gluon-fusion, H^-fusion, and btxf) contributions are summed) and 
decay in the 77 channel, for the same parameter choices as in fig. 15: 
a) h, LHC; b) H, LHC; c) A, LHC; d) h, SSC; e) H, SSC; f) A, SSC. 
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For the sake of comparison, the SM values are also indicated. QCD 
corrections to the gluon-fusion mechanism are not included. 

Fig. 21: Significance of the inclusive — > 77 signal, in the plane defined by m^ 
and a ■ BR(<p — > 77), for the CMS detector proposal at the LHC, with 

an energy resolution AE/E = [2%/ JE(GeV)] + 0.5%. The solid lines 

are contours of constant S/ yB, where S is the signal and B is the 
background. The dashed line corresponds to the SM Higgs, including 
QCD corrections to the gluon-fusion mechanism. Courtesy of C. Seez 

Fig. 22: Contours of constant cross-sections times branching ratios, in the (m^, tan /3) 
plane, for the inclusive — ► 77 channel: a) h, LHC; b) H, LHC; c) h, 
SSC; d) H, SSC. The choice of m t and rriq is the same as in fig. 15, and 
QCD corrections to the gluon-fusion mechanism are included. 

Fig. 23: Cross-sections for associated W<f> and tt<f> production, times branching 
ratios for the — > 77 channel, for the same parameter choices as in 
fig. 15: a) Wh, LHC; b) WH, LHC; c) Wh, SSC; d) WH, SSC; e) tth, 
LHC; f) ttH, LHC; g) tiA, LHC; h) tth, SSC; i) ttH, SSC; j) ttA, SSC. 
For the sake of comparison, the SM values are also indicated. 

Fig.24: Contours of constant L = [2a{ti(j)) ■ BR(t -> Wb) + a(W</>)] ■ BR(</> -> 
77) • BR(W — > Iv), for the same choice of mt and m^ as in fig. 15: a) h, 
LHC; b) H, LHC; c) h, SSC; d) H, SSC. 

Fig. 25: Cross-sections for inclusive H production (the gluon-fusion, P^-fusion 
and bhcf) contributions are summed) and decay in the Z*Z* — > 41^ 
channel (I = e, /x), for the same parameter choices as in fig. 15: a) LHC; 
b) SSC. For the sake of comparison, the SM values are also indicated. 
QCD corrections to the gluon-fusion mechanism are not included. 

Fig. 26: Contours of constant cross-sections times branching ratios for H — ► 
Z*Z* — > 4^, for the same choice of m^ as in fig. 15 and: a) m t = 
140 GeV, LHC; b) m t = 140 GeV, SSC; c) m t = 120 GeV, LHC; 
d) mt = 120 GeV, SSC; e) m t = 160 GeV, LHC; f) m t = 160 GeV, 
SSC. QCD corrections to the gluon-fusion mechanism are included. 
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Fig. 27: Cross-sections times branching ratios for — > r + r , for the same pa- 
rameter choices as in fig. 15: a) h, LHC; b) H, LHC; c) A, LHC; d) h, 
SSC; e) H, SSC; f) A, SSC. For the sake of comparison, the SM values 
are also indicated. 

Fig. 28: Contours of constant cross-sections times branching ratios for — > 
t + t~, for the same choice of m t and rriq as in fig. 15: a) h, LHC; b) H, 
LHC; c) A, LHC; d) /i, SSC; e) /J, SSC; f) A, SSC. The vertical lines 
in c) and f) correspond to rriA = 60 GeV. 

Fig.29: Contours of constant AR = [BR(t -> H + b)-BR(H + -> r+z/ r )]/[ J B J R(t - 
• J Bi?(H/+ -> for m 9 - = 1 TeV and: a) m t = 140 GeV; 

b) m t = 120 GeV; c) m t = 160 GeV. 

Fig. 30: Pictorial summary of the discovery potential of large hadron colliders 
for m~ q = 1 TeV and m t = 140 GeV. 
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